Magnetic excitations in NpCoGa 5 
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We report the results of inelastic neutron scattering experiments on NpCoGas, an isostructural 
analogue of the PuCoGas superconductor. Two energy scales characterize the magnetic response 
in the antiferromagnetic phase. One is related to a non-dispersive excitation between two crystal 
field levels. The other at lower energies corresponds to dispersive fluctuations emanating from the 
magnetic zone center. The fluctuations persist in the paramagnetic phase also, although weaker in 
intensity. This supports the possibility that magnetic fluctuations are present in PuCoGas, where 
unconventional d-wave superconductivity is achieved in the absence of magnetic order. 

PACS numbers: 71.27.+a; 75.30.-m; 78.70.Nx; 75.30.Et 
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Many complex and diverse phenomena have been ob- 
served in the series of isostructural RTX5 intermetallic 
compounds containing cerium or a light actinide element 
R, a transition metal T (Fe, Co, Ni, Rh, Ir), and an 
element of the boron group X (Ga or In). Current re- 
search on this family is enhanced by the possibility to 
explore links between unconventional superconductivity, 
magnetism, and non-Fermi liquid behavior in proximity 
of a quantum critical point [l| . Especially the discovery 
of d-wave superconductivity in PuCoGas at the surpris- 
ingly high critical temperature of T sc = 18 K \2\ raised 
much interest. No long-range magnetic order is observed 
in PuCoGas, but theoretical work [!, 0, [ft? ] as well as 
recent NMR experiments suggest the presence of spin 
polarization and hence strong magnetic correlations. Ev- 
idence for non s-wave superconductivity is reported also 
for CeRhlns and CeCoIns, where heavy fermion charac- 
ter is consistent with proximity to magnetic order []]]. 

Although magnetic fluctuations have been proposed 
to stabilize superconductivity for both Pu and Ce-based 
compounds, no direct measurement of magnetic fluctua- 
tion spectra has yet been reported. This is mainly be- 
cause of the large neutron absorption cross-section of 
239 Pu and 113 In, and the lack of large single crystals of 
242 PuCoGa5. Thus, at present, there is no information 
on the electronic excitations in this interesting set of ma- 
terials. 

Here we present the results of inelastic neutron scatter- 
ing (INS) experiments on the neptunium-based analogue 
NpCoGas. This compound (space group P4/mmm, a 
— 4.277 A, c = 6.787 A) shows antiferromagnetic (AF) 
order below Tjy = 47 K, with an ordered Np moment 
H = 0.84/is Q pointing along the tetragonal c-axis 
The magnetic structure is defined by the propagation vec- 
tor (0 1/2). This is exactly the reciprocal space position 
where theory predicts the dominant weight for the mag- 
netic fluctuations in PuCoGas ®- We therefore expect 
some similarities between the latter and the magnetic ex- 
citations in the antiferromagnetic compound NpCoGas. 



For this experiment we used a large single crystal with 
a mass of 1.1 g, grown and characterized at the Institute 
for Transuranium Elements, Karlsruhe. Measurements 
were performed with the thermal triple axis spectrome- 
ter IN8 at the Institut Laue-Langevin, Grenoble, using 
a Si(l 1 1) focussing monochromator, open collimation, 
two graphite filter for second-order contamination sup- 
pression, and a pyrolytic graphite (0 2) analyzer. 

The final wave- vector was fixed at kf = 2.662 A" 1 . 
The crystal was oriented with the a- and c-axis in the 
scattering plane. Several structural Bragg peaks were 
measured to test its quality, found to be adequate for 
an INS experiment. The rocking curve (0 sample rota- 
tion) at the (0 3) Bragg reflection shows a double-peak 
structure, the two maxima being separated by about 1.6 
degrees and each peak having 0.9 degrees full width at 
half maximum (FWHM). The sample orientation was op- 
timized for the crystallite with largest volume. Longitu- 
dinal scans at the (2 0) and (0 3) Bragg peaks gave 
Aq h = 0.026 and Aq t = 0.045 r.l.u. FWHM, respectively. 

Before searching for the magnetic inelastic response, 
we measured the neutron groups corresponding to the 
transverse acoustic (TA) phonons propagating along 
[0 1] to determine the lowest possible vibrational ex- 
citation energy at the magnetic zone center (h k 1 + 1/2). 
Constant-Q scans were performed around the (2 1/2) 
reciprocal lattice position with the sample kept at T = 
10 K. The results obtained are shown in Fig. [T] The 
observed energies and intensities are in good agreement 
with the results reported for UCoGas [10(. In particu- 
lar, the intensity at (2 0.7) is very weak in all cases. 
Since the TA[0 1] phonon branch reaches an energy of 
7.2 meV at the magnetic zone center, any excitation with 
energy smaller than this value cannot be of phononic ori- 
gin. 

Having established the energy of the lowest vibrational 
excitations, we looked for the magnetic response at the 
(1 1/2) reciprocal lattice point. At this position, the 
intensity of the TA phonon is expected to be weaker than 
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FIG. 1: Neutron groups observed with the sample at 10 K in 
constant-Q measurements at different wave vectors (2 q z ). 
The observed excitations correspond to TA phonons propa- 
gating along the [0 1] direction. 

at (2 1/2) by almost an order of magnitude, because 
of the geometric factor appearing in the neutron cross- 
section. The results obtained are shown in Fig. [51 Data 
were taken as constant-Q scans at different Q = (1 q z ) 
positions, and fitted to two Gaussian line-shapes. One 
Gaussian curve represents the contribution of the TA 
phonon; it is centered at the energy shown in Fig. [1] for 
the corresponding peak. The intensity scales with the 
square of the neutron momentum transfer, as expected. 
The second Gaussian curve is attributed to a propagating 
magnetic excitation, with an energy gap of 5.5±0.2 meV. 
Moving away from the magnetic zone center (1 1/2), 
the energy of phonons decreases whereas the energy of 
the magnetic peak increases. 

In addition, a weaker non-dispersive signal is observed 
at about 10.5 meV, which we attribute to a mean field 
excitation between crystal field (CF) levels. Due to its 
non dispersive behavior, this excitation is also visible in 
the scan performed at a generic reciprocal lattice position 
(0.8 1.2), as shown in Fig. H 

Equivalent scans at (2 1/2) and (3 1/2) establish 
an increase in the phonon contribution only, and confirm 
the magnetic origin of both, the propagating and the 
local excitation. 

To better characterize the non-local magnetic re- 
sponse, we followed its temperature dependence by per- 
forming constant-Q scans at the (0 3/2) reciprocal lat- 
tice point (Fig. [3]). At this position, the TA phonon con- 
tribution is minimized and the longitudinal acoustic (LA) 
phonon is expected at higher energies only. In addition, 
at such a small momentum transfer, all phononic signals 
are weak. The scan at the AFM zone center is compared 
with a scan at Q = (0.2 1.47), which was reached by a 
= 12.5 degrees sample rotation away from the optimal 
value for the (0 3/2) position. This position might be 
considered to reflect the Q-independent background (see 
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FIG. 2: Neutron groups corresponding to excitations propa- 
gating along [0 1] at 10 K. The solid lines are fits to the sum 
of two Gaussians and a sloping background, measured at Q 
= (0.8 1.2). The third peak at 10.5 meV (dashed vertical 
line) is attributed to a local, non dispersive excitation. The 
vertical tags indicate the energy of the TA phonons, and are 
labelled by the corresponding q z value. 

discussion below). By subtracting the background from 
the intensity at the AFM zone center and dividing by 
the Bose factor, we obtain the temperature dependence 
of the imaginary component of the dynamical suscepti- 
bility (Fig. [3]b). A phonon contribution is visible at 
8.4 ±0.6 meV and, as expected after Bose factor normal- 
ization, its intensity is temperature independent. On the 
other hand the magnetic intensity at 5.5 meV becomes 
weaker as the temperature increases (inset Fig (3]), but it 
is still clearly visible above the Neel temperature, Tn- 
The large difference between the intensities measured at 
(0 1.5) and (0.2 1.47) confirms that this signal cannot 
be attributed to a crystal field excitation, neither below 
nor above the Neel temperature Tjv. Indeed, the inten- 
sity of a crystal field excitation would change only as the 
square of the magnetic form factor, which is essentially 
the same at the two positions. 

It is interesting to compare the static and dynamic 
magnetic intensities at both the (1 1/2) and the 
(0 3/2) positions. The static magnetic moment points 
along the tetragonal c-axis, and therefore all Bragg reflec- 
tions of (0 L)-type have zero neutron intensity. Since 
the inelastic magnetic intensity at 5.5 meV of both reflec- 
tions is comparable, we suspect the dynamic fluctuations 
to show an (at least) significant transverse component. 

FigH] shows longitudinal Q-scans around Q = (0 q z ) 
and transverse Q-scans around Q = (q x 3/2), with 
neutron energy-transfer fixed at 5.5 meV. The non-local 
character of the 5.5 meV excitation at T = 10 K is 
confirmed by the strong Q dependence of the intensity, 
which peaks at the magnetic zone center. The peaks at 
(0 2.85) and (very much weaker) at (0 1.15) are due 
to the nearby LA phonon. At this stage of our investi- 
gations, with a relaxed instrumental momentum resolu- 
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FIG. 3: (a) Constant-Q scans at different temperatures at the 
Q = (0 3/2). A constant vertical offset of 1000 counts is ap- 
plied for clarity. The solid lines are scans at Q = (0.2 1.47). 
(b) Imaginary component of the dynamical magnetic suscep- 
tibility: squares 10 K, open circles 42 K, up triangles 52 K, 
down open triangles 100 K. The solid lines are fit to two Gaus- 
sian line-shapes centered at 5.5 ± 0.2 and 8.4 it 0.6 meV 
(vertical dashed lines), with intensities shown in the inset. 




FIG. 4: Q scans at constant energy transfer, hw = 5.5 meV. 
Left panel: longitudinal scan along the [001] direction, Q = 
(0 q z )- Right panel: transverse scan around (q x 1.5), 
performed by O rotation of the sample. Neutron intensity is 
plotted as a function of q x . The arrow indicates the position 
corresponding to the scan at Q = (0.2 1.47) in Fig. [3] A 
vertical offset of 0.2 counts/monitor is applied for clarity to 
data at different temperature. 



the excitation modes observed at 10 K in this experi- 
ment. The TA phonon dispersion falls in between the 
dispersion curves reported for the Pu and U homologous 
compounds [l(], 12 1- Also the point measured for the 
LA phonon is in good agreement with previous studies of 
the vibrational response in UCoGas and PuCoGas. The 
dispersion of the low-energy magnetic mode is typical of 
a spin-wave excitation in an antiferromagnet. The non- 
dispersive local excitation at 10.5 meV is also shown. 

Previous studies [12] assume the 5f electrons to be es- 
sentially localized and predict a crystal field (CF) level 
at about 8 meV in the paramagnetic phase, directly cou- 
pled with a doublet ground state. On this basis, an Ising 
mean field (MF) model has been proposed to describe 
the magnetically ordered phase [13j. The Hamiltonian 
can be written as the sum of a CF term, 



tion and without polarized neutrons, we cannot establish 
the origin of the large background signal, which increases 
with temperature. This prevents any quantitative anal- 
ysis of the peak width and therefore the spatial coher- 
ence length of the dispersive magnetic excitation. Wc 
note, however, that at 10 K the magnetic peak has a 
width comparable to that of the LA phonon, suggesting 
a long-range character. Approaching Tjv, the magnetic 
response becomes much broader in Q, and only a short- 
range in-plane correlation appear to survive. Polarization 
analysis experiments are planned to address this point. 

Fig. [5] shows the dispersion relation along [0 1] of 



H CF = BlOl + BlOt + B°Og + BfOi + B%0% (1) 

a nearest-neighbors exchange interaction within the a- 
b planes, 



<ij> 

and an interlayer exchange term 
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FIG. 5: Dispersion relation along [0 1] for the TA phonon 
(full circles), the LA phonon (open circle), and the two mag- 
netic modes (full and open diamonds). Data taken at T = 
10 K. Dashed and dotted lines are phonon dispersion curves 
reported for PuCoGas and UCoGas, respectively [ToL Hl|. 

where 0\ are Stevens operator equivalents, B\. the CF 
parameters, A 2 and Ajr are the intralayer and interlayer 
exchange constants, respectively, i,j label the Np site in 
the a-b plane, and I indexes the layers. 

Following Aoki et al. [12|, the lowest lying CF eigen- 
states have the form 

\d±) = a|±3> + Vl-a 2 |Tl> 
\s) = c(|+4) + | -4)) + v / l-2c 2 |0>. (4) 

The doublet \d±) and the singlet \s) are separated by 
a gap Acf, whilst the remaining states of the 5 14 multi- 
plet have much higher energies and do not affect the low 
temperature properties. The exchange Hamiltonian is 
treated in MF approximation and the self-consistent so- 
lution is sought. The coefficient a is fixed by the value of 
the ordered moment, and the sum A z + A^ is varied in or- 
der to fix Tj\r to the experimental value. A non-dispersive 
INS transition at 10.5 meV is indeed allowed if Acf = 
6 meV. This bare CF gap is close to the value proposed 
in [12( 1 to explain the bulk properties in the paramag- 
netic phase. In the ordered phase, the doublet would 
be split by A ex = 9 meV, however the corresponding 
INS transition would not be dipole-allowed. Therefore, 
within this model, the dispersive excitation appearing 
at 5.5 meV at the magnetic zone center cannot be un- 
derstood. Its observation indicates that a fully localized 
dipolar model fails to describe the dynamical magnetic 
response of NpCoGas . 

In conclusion, our inelastic neutron scattering exper- 
iment gives evidence for the existence of two energy 
scales for the magnetic excitations in the ordered phase 
of NpCoGas. One is related to a non-dispersive exci- 
tation between two crystal field levels (with a bare gap 



modified by mean- field exchange), and a lower one corre- 
sponds to dispersive spin fluctuations. The CF splitting 
is compatible with earlier predictions based on essentially 
localized dipolar models, which reproduce the bulk mag- 
netic properties in the paramagnetic phase [12| . The dis- 
persive magnetic response, on the other hand, reveals the 
presence of non-local excitations with an energy gap very 
close to that deduced from the low temperature depen- 
dence of the specific heat [8|. 

As shown in Fig. [3l an inelastic magnetic response is 
visible also at 52 and 100 K. This signal cannot be at- 
tributed to a CF transition, because of the different inten- 
sities measured at Q = (0 3/2) and at Q = (0.2 1.47). 
Therefore, although broader and weaker, the low energy 
excitation resulting from the spin-spin temporal correla- 
tion clearly persist in the paramagnetic phase. This sup- 
ports the possibility that magnetic fluctuations centered 
around (0 1/2) are present, as suggested by theory [f|, 
also in PuCoGas, where they have been proposed as a 
possible mediating boson responsible for the appearance 
of superconductivity [l4T | . 

Finally, it appears that the 5/ electrons in NpCoGas 
are neither completely localized nor itinerant. This dual 
nature is typical of the 5/ shell in intermetallic com- 
pounds and its elucidation is one of the main current top- 
ics in the physics of strongly correlated electron systems. 
Recent Mossbauer spectroscopy studies of NpTGas point 
to a tendency towards localization of the 5/ states with 
increasing number of d electrons on the transition metal 
T [15J . It will therefore be interesting to extend our INS 
investigation to NpFeGas, where an essentially itinerant 
response is expected, and to NpNiGa 5 , where the local- 
ized dynamics should dominate. 
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